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Four-quark condensates and chiral symmetry restoration in a resonance gas model 

Stefan Leupold 
Institut fiir Theoretische Physik, Universitat Giessen, Germany 

As an alternative to the two-quark condensate we propose a specific four-quark condensate as 
an order parameter of chiral symmetry restoration. We show that this four-quark condensate is 
closer connected to observable quantities. Within a resonance gas model we calculate the in-medium 
changes of two- and four-quark condensate as functions of temperature and baryo-chemical potential. 
In this way we estimate the line of chiral symmetry restoration in the temperature-potential plane 
and also as a function of energy and baryon density. It turns out that the line determined from the 
vanishing of the four-quark condensate is extremely constant as a function of the energy density. 

^- . I. INTRODUCTION 

o '■ _ 

^^ ' In QCD, the theory of the strong interaction, the masses of the lightest quarks are very light ||l| as compared to 

Cn ' typical hadronic scales. Therefore, QCD has an approximate chiral symmetry. From the absence of degenerate chiral 

?H partners and from the existence of very light pseudoscalar mesons one concludes that chiral symmetry is realized in 

^^[ the Nambu-Goldstone phase, i.e. it is spontaneously broken (see e.g. |2| and references therein). Consequently, the 

light pseudoscalar states are interpreted as Goldstone bosons. 

vQ ' It is expected that in a hot and dense strongly interacting medium chiral symmetry gets restored. To quantify the 

degree of restoration of a symmetry one typically uses an order parameter. Indeed, concerning chiral symmetry several 

7—i ■ quantities are conceivable as an order parameter. There are at least two considerations vi^hich decide how useful a 

^ [ possible order parameter actually might be, namely a) observability and b) simple connection to QCD. Clearly, it 

OO . would be preferable to have a quantity which was observable in vacuum and also in a strongly interacting medium. 

In this way one could experimentally test the process of chiral symmetry restoration. Unfortunately such a "smoking 

gun" is hard to find (cf. e.g. Q)- But it should be clear that it is appreciable as a first step, if the order parameter is 

f~^ ■ at least measurable in the vacuum. From the theoretical point of view it would be preferable to have a quantity which 

\^ fulfills criterion b), i.e. which can be easily expressed in terms of quark and gluon fields and which can be determined 

f^ in lattice QCD Q or quark model calculations. Basically this criterion b) boils down to the requirement that the order 

^^ ' parameter is a condensate. It is the purpose of the present paper to propose a specific four-quark condensate as a 

f^ ; very useful order parameter. First, however, we shall discuss order parameters which are more frequently considered. 

I ' We will especially comment how well criteria a) and b) are met. 

Obviously, the appearance of Goldstone bosons is closely connected to spontaneous symmetry breaking. Therefore, 
it is natural that one can find an order parameter of chiral symmetry restoration which is connected to the properties 
of Goldstone bosons. Indeed, the pion decay constant F^^ is such a quantity. It vanishes, if chiral symmetry is restored. 
Clearly, i^^ can be extracted experimentally from the life time of the pion. In principle, -FV is connected to a quark 
current, but the expression is non-local. In addition, it is not so easy to determine it from lattice QCD (cf. e.g. [3 
5^ ' and references therein). Thus, criterion a) is met by the pion decay constant, but criterion b) only in part. 
Cu I Frequently used order parameters are the two-quark condensates (uu), (dd) and (ss) or the corresponding two- or 

three-flavor averages.^ In lattice QCD calculations a dramatic rise of the pressure as a function of the temperature 
has been observed at a specific temperature Tc- This points towards a phase transition or at least (and more likely) 
a rapid crossover. It has been found that the two-quark condensate also shows a drastic change at the very same 
temperature Tc (more precisely, the corresponding susceptibility peaks at this temperature) |6|. This suggests that 
the two-quark condensate is a proper tool to study chiral symmetry restoration. On the other hand, it might be that 
the situation is different for non-vanishing baryo-chemical potential, i.e. at finite baryon density. Indeed, one can 
imagine scenarios where the two-quark condensate vanishes while chiral symmetry is still broken (see e.g. Q)- On top 
of that problem, criterion a) is not met by the two-quark condensate: It is not directly connected to an observable 
quantity. Here a closer look is useful: In the standard scenario of chiral symmetry breaking (xSB) one gets Q (see 
also 13) 

F^M^ « - {rriu + md){uu) « -(m„ -f md){dd) (1) 
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^ For the following qualitative discussion we will simply talk about the two-quark condensate without specifying the flavor. When it comes 
to quantitative statements we will specify which two-quark condensate is considered. 



(where isospin-breaking effects are neglected for the condensates). We see that it is not simply the two-quark con- 
densate but its combination with the quark mass which is connected to observables^ — and also here not to one 
observable, but to a combination of two of them, the pion decay constant and the pion mass M^^. From Q it is 
obvious that it is possible that the two-quark condensate vanishes with a pion decay constant which stays more or less 
constant. A change of the pion mass would be sufficient to change the condensate. In such a case, chiral symmetry 
would still be broken. Besides (^ there is so far no relation derived from first principles which directly connects 
two-quark condensates with observable quantities. 

To summarize, the pion decay constant satisfies criterion a), but no so much b), whereas it is just the other way 
round for the two-quark condensate. Therefore, it is reasonable to explore also other order parameters of chiral 
symmetry breaking. Note that not much lattice QCD studies have been performed here. The mixed quark-gluon 
condensate has been studied on the lattice in [IJ. We will briefly comment on this condensate below f section llV)l . In 
the present paper we will mainly concentrate on a specific four-quark condensate. In the next section we will motivate 
in detail why this condensate qualifies as a useful order parameter. So far there are no lattice QCD calculations 
available for the temperature dependence of this four-quark condensate. In addition, an extension of lattice QCD 
to finite baryo-chemical potentials is still complicated (cf. e.g. [l3i[l3)- Therefore, we will evaluate the in-medium 
expectation value of the four- and also of the two-quark condensate within a resonance gas approximation in sections 
Illll and llVI The numerical results are presented in section A summary and an outlook are provided in section IVTI 
In an appendix some formal issues are discussed. 

II. GENERALIZED WEINBERG SUM RULES 

The two main arguments which point towards spontaneous breaking of chiral symmetry are the appearance of light 
states (interpreted as Goldstone bosons) and the absence of degenerate chiral partners. In the following we focus on 
the second aspect. 

We start with the retarded {R) current-current correlator [lj| 



U,Aq) ^ijd'x e^^^Rj; {x)jI (0)) (2) 

where j)f is the electromagnetic current with the quantum numbers of the p-meson, 

il = 2 (""^m" - dl^,d) . (3) 

The expectation value for an arbitrary operator O is defined by 

Tr(e>e-'3(^-'^^)) 
^^) ■" Xr(e-'5(//-pAr)) ('^) 

with the Hamiltonian J?, the baryon number operator N , the inverse temperature /3 = 1/7^, the baryo-chemical 
potential jjl and the trace Tr is taken over all possible n-body states. 

In a first step, we will decompose the Lorentz tensor n^j^(g). In the vacuum we only had the structures g^^ and 
g^gj, at hand. An equilibrated medium with finite temperature and baryo-chemical potential introduces an additional 
Lorentz vector n^ which is conveniently chosen to be proportional to the (conserved) baryonic current and normalized 
to 1. In the rest system of the medium it is simply given by n^ — (1,0,0,0). Since the current ^ is conserved by 
the QCD equations of motion we have 

g^n^,(g)-0. (5) 

We can construct two independent projectors Lfj_,y{q) and T^^{q) which both still satisfy current conservation 
g^L^,(9) = q^T^M = (cf. e.g. [H): 

r I . _ ((g ■n)qt,~ q^ n^) ((g ■ n) q^, - q^ n^) 

^M-(9) - g2 (q2 _ (^ . „)2) ' i^^) 

Ti,u{q) = gt^u- -^ - L^,^{q). (6b) 



Indeed, only this combination is renormalization group invariant lid . 



Hence, the correlator can be decomposed in the following way: 

n^, = n^T^, + n^i^, . (7) 

The scalar functions 11^ and H^ in general depend on q^ and (g • nf' . In the rest frame of the medium one can 
equally well state that 11^'^ depends on q^ and q^ — {q-n)'^ —q^ with the three- momentum g relative to the medium. 
The projectors T^;^ and I/^^ introduced in IJBJ are chosen such that they project on states which are transverse and 
longitudinal, respectively, with respect to q. 

For simplicity we restrict ourselves in the following to the case of vanishing (f, i.e. where the current is at rest 
relative to the medium. In this case there is no distinction between longitudinal and transverse states, i.e. 11-^ and 
n^ become equal: 

n^iq^q' - 0) = ^^(<Z^g2 ^ 0) ^: -q^R'^iq') . (8) 

R^ has a direct physical meaning in the time-like region s — q^ > 0. It is related to the cross section e'^e^ — > 
hadrons with isospin 1 via lOj 



127rlmi?*' . (9) 



a^ ^(e+e ^hadrons) 
a{e+e^ -^ ^+^^) 

In addition, it is related to the decay of the t into a neutrino and an even number of pions |l6tll7l|. 

If chiral symmetry is manifestly realized, i.e. not spontaneously broken, the current-current correlator for the vector 
current ||2Jl and the corresponding one for the axial-vector current 

j;f = 2 (■"^A'75U - ^7^75^) (10) 

must be degenerate. In other words, in this case the spectrum of hadronic states which couples to j^ must be 
degenerate with the corresponding spectrum which couples to j . These are the mentioned chiral partners. If chiral 
symmetry is spontaneously broken, there are no degenerate chiral partners and the following Weinberg sum rules |lq | 
hold (in the chiral limit) 

- fds {lmR^{s)~ImR'^{s)) = F^ , (11a) 



oo 

- fdss {lmR^{s)-ImR^{s)) = 0. (lib) 



The quantity R^ is connected to the vector correlator according to |(SJ). R-^ is the corresponding quantity for the 
axial-vector current IjlUI) with one subtle but important difference: In principle, also the pion contributes to R with 
a strength determined by the pion decay constant Ft^. This contribution is taken out from R and displayed explicitly 
in Hlla|l . Actually both quantities R^ and R^ can be observed in r decays, if one triggers on an even (for R^) or 
odd (for i?'^) number of pions in the final state p!a.[l7l|. 

The sum rules (|ll|l have been derived prior to the advent of QCD by using current algebra [lal- In other words, 
they are independent of QCD. In principle, one could also weight the difference R^ — R^ with higher powers of s. In 
contrast to (|ll|l such integrals are specific to QCD since quark fields and the strong coupling are involved (see below). 

From a practical point of view the predictions l|ll|) cannot be tested as they stand, since the integrals cover the 
whole range up to infinity. On the other hand, we expect that for large s the inclusive cross section in (Q and therefore 
R^ should be given by perturbative QCD [1^. Since perturbative QCD does not show chiral symmetry breaking we 
expect the same behavior for the axial-vector correlator, i.e. 19, 20] 

Imi?^(s) = Imi?^(s) = — fl + — ) for s > Sq (12) 

with a threshold value sq- 



Using (|12(l we can restrict the integration limits in (|ll|l and get 

So 

- fds {lmR^{s)^lmR^{s)) ^ F^ , (13a) 



So 

- jdss (lmi?^(s) -Imi?^(s)) = 0, (13b) 



So 

^ fdss^ {lmR''{s)-ImR^is)) - -Ina^iO^sB) ■ (13c) 



We have supp lemented the Weinberg sum rules by a third sum rule H13c|l which — as already announced — is specific 
for QCD '2(j\ . It involves the four-quark condensate 

(O^sb) = ((fl7,.75A'^u-(J7,.75A'^d)'-(w7^A''w-J7^A"d)2> (14) 

where Xa denotes the color Gell-Mann matrices. For an in-medium generalization of ( 1 1311 see also |2l] |. 

Actually the quantity as(C'xSB) is not a renormalization group invariant object ^T\. A serious derivation of (|13c|) 
reveals that it must be evaluated at the scale Sq (cf. [23 and references therein). There are also radiative corrections 
which we have not displayed explicitly. Since the time when four-quark condensates were introduced 12 0| it is a 
lively discussed issue whether they can be related to the two-quark condensate by factorization (see e.g. J2J| and 
references therein). This discussion is complicated by the fact that four-quark condensates are difficult to determine 
experimentally. In addition, also the size of the two-quark condensate is not so well determined: With (Q one can 
relate the two-quark condensate to the very well known pion mass and decay constant and to the much less known 
current quark masses [ij. This uncertainty doubles since under the assumption of factorization the square of the 
two-quark condensate is connected to the four-quark condensates. With all these problems in mind we note that it 
has been found recently in 25j that in vacuum the four-quark condensate ()14|l indeed factorizes. This finding is based 
on the recent data on r decays and on a technique which we will describe in a moment. On the theory side there is 
one particular limit in which vacuum four-quark condensates are connected to the two-quark condensate, namely for 
a large number of colors Nc J2g . For in- medium condensates this is again more subtle 24] . We will come back to 
that point below. 

From a phenomenological point of view the sum rules H13|l are still unsatisfying: As already mentioned, nowadays 
experiments can address the difference of correlators in t decays. Therefore, one can explore the integrals of H13|l 
only up to energies below the t mass. It turns out that the left hand sides of the sum rules are not convergent as a 
function of the threshold sq for values below the t mass (squared). By proper combinations of (|13|) one can, however, 
decrease the sensitivity on sq and obtain sum rules which show a much better convergence behavior already at the 
experimentally accessible energies [23, 123 : 



so 

— [ds(so-s) (Imi?^(s)-Imi?^(s)) = F^ , (15a) 
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So 

^ fdssiso - s) (Imi?^(s) - ImR^is)) = i 7ra,(OxSB) • (15b) 



One can obtain the left hand sides of these sum rules from the measured spectra determined from r decays. The 
right hand side of (|15a|l is given by the experimentally measured pion decay constant. It has been found ,2a| that 
both sides of H15a() agree very well. In the same way, the left hand side of (|15b|) has been determined. Here it has 
been found that the result agrees well with the conjecture that the four-quark condensate H14|l factorizes into the 
square of the two-quark condensate |23. This is exactly what large- TVj, considerations suggest. Below we will take the 
finding of !25l | as a motivation to explore also the in- medium changes of the four-quark condensate (|14|l in the large- TVc 
approximation. Note that we do not claim here that an arbitrary four-quark condensate is always connected to the 
square of the two-quark condensate. First of all, as we will see below, even for the fo ur-q uark condensate at hand, 
its in- medium change is not given by the change of the two-quark condensate (see also J24|). Second, from the point 
of view of chiral symmetry breaking, it makes sense to connect a chirally odd four-quark condensate, e.g. (|14|l . to the 



two-quark condensate. Concerning chirally invariant four-quark condensates^, on the other hand, it is not so clear 
whether there should be any deeper connection to the two-quark condensate. In the present work we concentrate on 
one specific four-quark condensate, namely (|14|l . In vacuum, this condensate seems to factorize [25j . 

From (|15|l it is obvious that both quantities F^ and the four-quark condensate (|14|l can be regarded as order 
parameters of chiral symmetry breaking (xSB). They must vanish if chiral symmetry is restored. To be precise, in a 
medium the relevant definition of Fj^ is the coupling of the pion to the temporal component of the axial- vector current 
(|10|l 28]. Concerning their role as order parameters, we want to stress again that F^ and as(0^sB} are actually more 
advantageous as compared to the two-quark condensate (qq) which is the "standard" order parameter: The former 
are connected to quantities which are in principle measurable, namely the left hand sides of relations (|15|l .^ 

The in-medium changes of F^ are addressed in J23,|30| for a non-interacting pion gas (density p^) and e.g. in [23 
for a non-interacting (and cold) Fermi sphere of nucleons (density Pat). For finite temperatures (pion gas) one gets 



F^iT)=F^^{l-^^ 



with the (scalar) pion density 



Ptt = 3 



d^k 



(27r)3 2Ek 6^"/^ - 1 



M^^O 1 j,2 

' 8 



(16) 



(17) 



and the pion energy Ek = yk'^ + M^. For cold nuclear matter (nuclcon Fermi sphere) the change of F^ can be 
related to some low-energy constants of the chiral pion-nucleon Lagrangian. For details we refer to |23|. Numerically 
the result is 



FUpn) = F^ 1 



Pn 



(0.52 ±0.08) 



(18) 



where po denotes nuclear saturation density. These results can be contrasted with the corresponding changes of the 
two-quark condensate {uu + dd) |29l l3ll | : 



{uu + dd) 



pionic med. 



{uu -\- dd)^ 



1- 
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The nucleon sigma term 



{uu + dd) 



nucl. inccl 
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Pn<^n 
F^ml 



1 Pn 
3 Po 



dm 



N 



dnin 



45 MeV 



(19) 



(20) 



(21) 



which appears in (|20|l can be obtained within the framework of chiral perturbation theory from low-energy -kN 
scattering ^33]. In principle, there are corrections to these linear-density results (|16|) . (|18|l . H19() and (|20|) which are of 
higher powers in the respective density. They involve correlations between the constituents of the medium. For the 
quark condensate in a pion gas this has been worked out in |3lj| . For all other cases model independent predictions 
are hard to make for the terms beyond linear order in the density. 

Comparing (|16|l with (|19|l or (|18|l with (|20|l shows that a bold extrapolation of these formulae to the point where 
the respective order parameter vanishes would lead to different critical densities. This indicates a breakdown of the 
linear-density approximation for at least one of the compared respective quantities, probably for both. 

Beyond these special cases of a (non-interacting) pion or nucleon gas the changes of F^ are hard to estimate (except 
for specific models). The situation is somewhat better for the four-quark condensate 1)14(1 to which we turn for the rest 
of this work. We will compare the obtained results to the in-medium changes of the two-quark condensate {uu + dd). 
We note in passing that the in- medium behavior of {ss) is different from the one of {uu -\- dd) J3J| . Since we started 



out from the non-strange correlators R and i? , a comparison to 



dd) is most appropriate. 



^ which appear e.g. in QCD sum rules for the u) meson |27|l 

* In practice this is a very tedious task: Access on BX can be obtained from dilepton spectra as discusse d above ; see also relation 1^ . 

Indeed, one tries to measure BX also in a medium (see e.g. Q). In vacuum, R is related to r decays |lfit llTl as already discussed. 

This is obviously difficult to measure in a medium. 



III. RESONANCE GAS MODEL 

The in- medium changes described above in equations (|16() . H18() . H19() and l|20() have been determined under the 
assumption that the respective medium is described by a non- interacting gas of the respective most abundant particles 
(pions for a medium with finite temperature and nucleons for cold nuclear matter) . Of course, to describe a medium 
with finite temperature and finite baryo-chemical potential one can easily generalize the previous results by using a 
gas of pions and nucleons (and anti- nucleons).^ Clearly, this is not the full story — except for low particle densities. 
Indeed, the in-medium expectation value Q is defined with respect to all possible states. Beyond the consideration of 
single-particle states one has to take into account n-particle states with n > 1, e.g. two-pion, pion-nucleon, nucleon- 
nucleon, three-pion states and so on. The influence of n-particle states is accompanied by higher powers in the 
respective densities and therefore suppressed for low temperatures and low baryo-chemical potentials. This justifies 
the previous results as the respective leading terms in a low-density expansion. However, we want to understand now 
how to go beyond this linear-density approximation — at least approximately. 

One might wonder why we only talked about pions and nucleons so far and not about their excitations. Indeed 
we will come to the excitations soon. From a principal point of view, however, we stress that there is no need to 
include them in the traces which appear in Q): For a complete set of states it is sufficient to consider all single- and 
many-body states built from the stable states (for flavor SU{2): pions, nucleons and anti-nucleons). As demonstrated 
long time ago in J35j the (unstable^) excitations emerge from the scattering phase shifts which come into play when 
considering many-body states in the traces in Q). E.g. a hadron resonance which is formed in pion-nucleon scattering 
is automatically taken into account when one considers the contributions of two-particle states (pion-nucleon) in the 
traces in (0J. 

On the other hand, it is in practice intractable to calculate the contributions of all n-body states to the expectation 
value Q. Now we can turn the previous argument around: Physically we expect that the contribution of a many- 
body state is the more important, the larger the many-body correlation is. Clearly, a large correlation is found, if a 
hadronic resonance is formed. Therefore one might approximate the complete set of (stable) many-body states by a 
sum over all one-body resonance states |3g. Below we will evaluate the in-medium two- and four-quark condensate 
in this resonance gas approximation (see also e.g. |3lLl37l l3q for successful applications of that idea). 

The following aspects should have become clear from the previous discussion: 

a) The resonance gas approximation takes into account part of the many-body correlations of the stable states which 
form the medium. In that way one goes beyond the linear-density approximation discussed above. 

b) The resonance gas approximation is not a systematic expansion in powers of the densities of (stable) medium 
constituents. One merely tries to take into account the most important parts of the higher density terms. 

c) The resonance gas approximation can only make sense, if the many-body correlations are governed by hadron 
resonances. This excludes the application to systems with low temperatures and large baryo-chemical potential. In 
such systems nucleon-nucleon correlations are the most important ones. 

In the resonance gas approximation we flnd 



(O)med. ~ (0|O|0) + J2 Px{X\0\X) (22) 

X 

where |0) denotes the vacuum state and O an arbitrary operator, px is the scalar density of states X: 

(Pk mx 
(27r)3 'E^ 



. (Pk mx 
P^ ^ / 7T^ 'W~ "^f/bKEx - sfj.) (23) 



with 

np/BiE) 



exp ( I ) ± 1 



(24) 



Ex = y'nT'x + ^^ ^^'^ s = 1 for baryons, s = — 1 for antibaryons and s = for mesons. The ± sign and the label 
F/B refer to baryons and mesons, respectively. We sum over all states X identified by the particle data group 0]. 



^ For our qualitative discussion we restrict ourselves to flavor SU{2). For our quantitative analysis we will include strangeness. 
® with respect to the strong interaction 



The normalization of the state \X) is chosen to be 



^x '^ ^3i 



{X(k)\X(k')) = —^(2nyS(k~k'). (25) 

mx 

We aheady note here that this norniahzation is appropriate for heavy states since it has a proper non-relativistic 
Umit. On the other hand, it is not quite adequate for Goldstone bosons where we would like to have states with a 
properly defined chiral limit. Of course, one can choose a different normalization, if one changes the definition of the 
density (|2^{|l accordingly. We will come back to that point below. 

We have to discuss the evaluation of H22|) separately for the two- and four-quark condensates. Concerning the type 
of hadron X we also will distinguish in the following between Goldstone bosons and other hadrons. If X is a Goldstone 
boson, the expectation values {X\uu + dd\X) and {X\0^sb\X) can be calculated using current algebra [l4ll39l|. We 
postpone this discussion to section Hvl 

To calculate the contribution of non-Goldstone states to the in-medium part of the two-quark condensate we follow 
the approach of jsj and generalize it to finite baryo-chemical potential (see also |13illl|)- We approximate^ 

{X\uu + dd\X) « {X\u^u + dU\X) = ll-ff' ^^'^°^'' (26) 

^ ' I / \ I ' ' ] 2 — Ns ioT mesons, ^ ' 

where Ng denotes the number of strange quarks in the state X. For the sigma term [S^l this amounts to the 
approximation 

_ J (3 - Ns) ruq for baryons, , . 

'^^ " \{2-Ns)mq for mesons. ^ '^ 

For the nucleon and the A(1232) we can get an idea how good this approximation is: Using niq « 7MeV and H21|l 
we find that 127|l underestimates a^ by about a factor of 2. On the other hand, a recent determination of cta yields 
20.6 MeV 42] in good agreement with (|27|l . We will explore the uncertainty induced by the estimate H26|l by an 
additional calculation where we simply increase ax by a factor of 2 for each particle species. The in-medium two- 
quark condensate (but not the four-quark condensate) has also been calculated within a resonance gas approximation 
in [i^. There an estimate different from (|26|l has been used for the quark condensate within a resonance. We will 
come back to this point below. 

Next we have to evaluate the expectation values appearing in (|22(l for the four-quark condensate H14|l . For the 
Goldstone bosons we note again that this can be performed using current algebra as will be discussed below in section 
IIVI For all other hadrons it is difficult to get an estimate for {X\OxSb\X). As already noted after equation (|14(l it 
is even a problem to pin down the vacuum expectation values of four-quark operators. For (0|O;^sb|0) it has been 
shown in |25j that it factorizes in the vacuum. In line with that finding, we use in the following the large- TV^ expansion 
p3Ll44l | where N^ denotes the number of colors. For the operator of interest, O^sBj we obtain |23.l45| 

(0|OxSb|0) = 0{N^) = 8 {0\uu\0)^ + o{N,) ; (28) 

{X\0^sb\X) =0{NJ) =%{Q\uu\{)){X\uu + dd\X)+o{Nc), (29) 

if X is a baryon; 

{X\0^sb\X) = o{N,) , (30) 

if X is a meson. 

Therefore, it seems that to leading order in \/Nc we only must consider baryonic states and we can relate the 
in-medium four-quark condensate to two-quark condensates. However, there is an additional implicit iVc-dependence 
in px which enters (|22|1 . We have shifted the discussion of that issue to appendix 1X1 The outcome of this discussion 
is that one should keep the respective leading order in 1/Nc for every type of hadron to get a serious estimate for all 
regions of T and ji (as long as T is not too small, cf. point c) in the discussion above). 

To obtain the respective leading 1/Nc contribution for a specific type of hadron we need a closer look at the 
derivation of (|29|l and (|30|l . We will summarize the essential points of [2J| and apply it to our case of interest. First 



Note that this relation becomes so simple using the normalization 1251 . 
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we note that the proof for (|29|l and H30() can be given for four-quark condensates involving white (cofor singlet) quark- 
antiquark operators. On the other hand, our condensate of interest (|14ll involves color octet operators. However, we 
can use a Fierz transformation to obtain 

{{u'J^'y^XaU - djfj^'ysXadf - {u-f^j,XaU - djf_,Xadf) = 

2 {{uu + dd) ) + 2 {{uij5U — dij^d) ) — 8 (uijsd di^^u) + 2 {{ui'f^u + di"f^d) ) + 2 {{uu — dd) ) — 8 { ud du) 




2 - 

- 1^ ((u7m75U - djf,-/5df - {ujf,u - d-y^df) . (31) 

JVc >. „ ' V ^ ' 

For an easier discussion and for orientation we have attributed meson states to the respective operators which have 
the same quantum numbers. The four-quark operators which appear on the right hand side of H31|) all have the generic 
structure qTq qT'q where T and V denote spinor and flavor, but no color structure. It has been shown in [2J| that for 
this type of four-quark operators the following holds for expectation values with respect to single-particle states X: 

{X\qTqqT'q\X) - (X|gTg|X) (0|qT'g|0) + (0|qTq|0) (X|qT'g|X) 

+ {X\qTq\0) {{)\qT'q\X) + {0\qVq\X) {X\qV' q\0) + terms subleading in l/N^. (32) 

Obviously, the first two terms on the right hand side of (|32l) can contribute for every type of X , but only if F and F' 
have the quantum numbers of the vacuum. This is only fulfilled by the first term appearing on the right hand side 
of (|31(l . We can use again (|26(l to evaluate the respective contribution. This is the factorized part of the four-quark 
condensate. 

Next we turn to the third and fourth term given on the right hand side of H32|) . These terms are contributions 
beyond factorization. First we note that these terms do not exist, if X is a baryon: A quark-antiquark operator qTq 
cannot create a baryon X from the vacuum. This simple fact finally leads to (|29|l . If X is a meson, the third and 
fourth term given on the right hand side of (|32|l only exist, if the quantum numbers of the meson match to the ones 
of qTq and qT'q. If this is the case, we need to know the overlap of the meson X with the quark currents qTq and 
qT'q. In principle, such a question can be addressed within lattice QCD. However, for higher excited hadrons this is 
a difficult task. In principle, we would need this information. In practice, however, we will find that we can safely 
neglect the respective contributions. To see how this comes about, let us first discuss which hadrons are of relevance 
here: In view of (|31|l we see that particles which yield contributions in leading 1/Nc order on top of the factorized 
part have quantum numbers of a) Goldstone bosons, b) /o or c) oq. We disregard the vector and axial- vector mesons 
here since their contributions are 1/Nc suppressed in (|31|l . Concerning the group a) we note that only the excited 
states are of concern here: Goldstone bosons will be fully taken into account in section Hvl There, no approximations, 
especially no large-A'c expansion is used besides the current algebra technique. Nonetheless, we note that we have 
checked explicitly that the current algebra calculations are not in contradiction to the large-A'c considerations used 
here. Still we have to be concerned with excitations of the Goldstone bosons and with /o and ao states. All these 
states have masses of about 1 GeV or higher |y.* As we will show below (cf. figure Q} these higher excited meson 
states are irrelevant from a practical point of view: For temperatures below, say, 200 MeV relevant contributions to 
the sum (|22|l over all resonances come from the low-lying mesons (pseudoscalar and vector nonet) and from baryons. 
Higher lying mesons are strongly Boltzmann suppressed and — in contrast to high lying baryons — do not have 
large degeneracy factors. It turns out that the sum of all baryon resonances — albeit also Boltzmann suppressed — 
adds up to a non-negligible contribution whereas the sum of all higher lying mesons remains irrelevant. Therefore, in 
practice we do not need the third and fourth term on the right hand side of (|32|l — except for Goldstone bosons to 
which we turn in the next section. Finally we note that in the spirit of the large-A'c expansion we also treat the rj' as 
a Goldstone boson (cf. e.g. (43 and references therein). 



We disr egard the /o{400 — 1200) which we consider as a loosely bound meson-meson molecule and not a genuine hadron resonance state 

mm. 



IV. GOLDSTONE BOSONS 

For Goldstone bosons we can use current algebra to calculate the contribution (X|0|X) in H22|l . First, however, 
we shall choose a normalization which is more appropriate for Goldstone bosons. Instead of H23|) and H25|) we use for 
Goldstone bosons |^3 

and 

(G(fc)IG(fc')) = 2^G {2TTf6{k-k') . (34) 

With this normalization (which has a well-defined chiral limit) we get using current algebra e.g. for pions |lj,|39j: 



1 

F3 



(^"|Ol7r^) = -^(0|[Q^,[Q^,O]]|0) (35) 



with 



i5 



i ld^xq{xhol5r''q{x) (36) 



and a Gell-Mann flavor matrix r°. 

Obviously this can be generalized to other Goldstone bosons as well. As already mentioned we treat the whole 
pseudoscalar nonet as Goldstone bosons in the spirit of the large- A^^, approximation. We use the decay constant 
corresponding to the type of particle in the denominator of (|35|) . For the flavor singlet we use t° = •\/2/3 1 . 

It turns out that t] and 77' do not contribute at all to (0)(SB)mod.- They do contribute to the two-quark condensate 
(see (I38|l and (|39l) below). Why there is no contribution to the four-quark condensate can be explicitly seen, if one 
does not use current algebra directly for (|14|l . but applies it separately to all the terms appearing on the right hand 
side of H31|l . i.e. after Fierz transformation. If A is an 77 or 77', the first operator on the right hand side of l|31|l 
contributes via the first two terms on the right hand side of (|32ll . The fourth operator on the right hand side of H31|) 
contributes via the third and fourth term on the right hand side of H32() . We have checked explicitly that in total 
these contributions indeed cancel. 

With relations like H35(l all expectation values can be traced back to vacuum expectation values. For the four-quark 
operator O^sb this leads to similar four-quark operators which in part deviate by their flavor content. We assume 
flavor symmetry of the vacuum to relate all resulting vacuum four-quark condensates to (OJOj^sbIO). No factorization 
(P5|l is needed here. 

With the same technique (|35|l the two-quark condensate {uu + dd)^cd. is related to the vacuum two-quark conden- 
sates {0\uu + dd\0) and (0|ss|0) (the latter appears using kaons instead of pions in pSfl ). Again we assume flavor 
symmetry of the vacuum: (0|ss|0} « {0\uu + dd\0)/2. 

Finally we want to comment on another order parameter of chiral symmetry breaking, the mixed quark-gluon 
condensate^ {qa^^G^'^q) with the gluon field strength G^" . In jllj it has been found that in lattice QCD the 
teinperature dependence of the two-quark and the mixed condensate are practically the same. It has been suggested 
in [ll| that this finding points towards a universal behavior of the order parameters. With the techniques presented 
here we can evaluate (g(T^i/G^''<z}med. for a gas of Goldstone bosons. It turns out that the operators Qf given in H3t)|) 
do not distinguish between 1, which appears in the two-quark condensate, and cr^vG^'^ , which appears in the mixed 
condensate. Therefore, for a gas of Goldstone bosons one gets 

(g-a^.G^'-g).„ed. ^ (0|ga^.G^'-g|0) ^^^^ 

{qq)mcd. {o\qq\o) 

On the other hand, for the four-quark condensate H14I) the temperature dependence is different (see below). Therefore, 
it is questionable to conclude only from the temperature dependence of the two-quark and the mixed condensate that 
there is a universal behavior. Note that we do not claim that a gas of Goldstone bosons is equivalent to a lattice 
QCD calculation. We just have presented an example where the temperature dependence of two-quark and mixed 
condensate agrees, but where other order parameters show a different behavior. 



For the following comment we do not distinguish between flavor SU{2) and flavor SU{3). 
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RESULTS 



Now we can take all terms together. We find for the four-quark condensate: 

CXD _ OO 

(CxSB)mcd. 



(OIOxSbIO) 





OO 

- T,2p2M2 Z^ (3 - N,)umB Jdk — npiEB - s^i) 

^ -^ B g 

OO 

/I 2 
dk — nsiEM) 



(38) 



where M denotes meson resonances except for tt, K , rj and 77'. The factor u parameterizes the uncertainty connected 
to the estimate (|26|l . In the figures below we will show calculations using u = 1 in agreement with (|26|l . In addition, 
we will use u = 2 to explore the uncertainty of our estimate (cf. the discussion after (|27|l '). 
For the two-quark condensate we get 



{UU + (id)nied. 

{0\uu + dd\0) 



oc oc 

3 f k'^ 1 f k^ 

dk -^ nsiE^) - ^_^^2 jdk — ns(^x) 



47r2F^2_y - p, 


00 



27r2F; 



K 



_L_ jdk |- ns{E,) ^ Jdk |- nniE,,) 

'0 

oc 

-J^Y.i^-^^)^^Bjdk^n,iEn^s^) 
B 

OC 

,^-^^^Y.i2~N,)umMJdk 



■nB{EM) ■ 



(39) 



For the numerics we use the parameter values listed in tabled 



quantity 


value in MeV 


ref. 


Fn 


92.4 


fll 


M^ 


140. 


fl] 


ruq 


7. 


m 


Fk 


113. 


\M 


Fr, 


124. 


[48] 


F,, 


107. 


[48] 



TABLE I: Parameter values used for actual calculations. 

As already noted we have treated the kaons as Goldstone bosons. Of course, one can study which formula we would 
have obtained, if we treated the kaon in the same way as the non-Goldstone bosons. This yields an opportunity to 
check our approximations. If kaons were considered as standard mesons instead of Goldstone bosons we observe first of 
all that there is indeed a relative factor 2 between the two kaon contributions in H38I) and (|39l) — this is what one would 
expect from factorization: In (|32|l the first two terms on the right hand side give the leading 1/A''c contribution, if X 
is a kaon. The third and fourth term drop out since the quantum numbers of the operators appearing in (|31l) do not 
match to the kaon. Of course, we see this same relative factor 2 in all non-Goldstone contributions in H39(l compared 
to (|38|l . Next we compare the kaon contribution in H39|l to the last contribution in H39|l : Instead of l/(27r^i^^) we 
would have for u — 1 (the factor 4 comes from the four types of kaons) 



2T:^FUa 



Mk. 



(40) 
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Indeed these two numbers deviate only by 6%. This finding gives us some confidence that our approximations are not 
too bad. 
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FIG. 1: Drop of the four-quark condensate 1381 (for u — 1) as a, function of the temperature T for three different baryo-chemical 
potentials: ^ = (full line), /i = 400 MeV (dashed), fi = 800 MeV (dotted). 



With (|38(l at hand we can study how the four-quark condensate changes in a medium. This is depicted in figure ^ 
We see that the four-quark condensate drops both with temperature T and baryo-chemical potential /i. Of course, for 
the two-quark condensate there emerges a plot qualitatively similar to the one shown in figure ^ We will not display 
these results, but concentrate in the following on the line where the respective condensate vanishes. This is depicted 
in figure 121 From inspecting H38|l and 1)39(1 one can already deduce that the four-quark condensate drops faster than 
the two-quark condensate and that the drop increases with the uncertainty parameter u. We observe these qualitative 
features also in figure |21 Quantitatively we deduce first of all that for vanishing chemical potential all lines yield a 
reasonable transition temperature between about 165 and 200 MeV. Without a bias which order parameter might be 
preferable we can view the lines in figure|21as estimates for the transition line to a chirally restored state. The highest 
and the lowest line mark the uncertainty of the estimates. On the other hand, as discussed in great detail in section 
^we prefer the use of the four-quark condensate as a quantity which is closer connected to observables. As we will 
see next, the transition line extracted from the vanishing of the four-quark condensate shows an additional interesting 
feature. 

An appealing aspect of the resonance gas approximation is the fact that one can easily calculate the energy and 
baryon density which correspond to given temperature and baryo-chemical potential: 



£ = 



E 



d3fc 



Eb npiEB - s^) + 22 



, (27r)3^ "'^ ^' ^J (27r)3 



(fk 



EMnsiEM) 



-^ ^ \dk k^Es [npiEB -fi)+ npiEB + m)] + ^ E M ^^^^f ^bI^m) 

B.noB h M n 



(41) 



and 



P = 



OC 

—2 ^ jdkk^ [np{EB - y) -nF{EB+ p)] 



(42) 
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FIG. 2: Vanishing of two- and four-quark condensates in temperature-potential plane for different values of u. Lines enumerated 
from top to bottom. Line 1: Vanishing of two-quark condensate (I39I I with u — 1. Line 2: Vanishing of two-quark condensate 
with u = 2. Line 3: Vanishing of four-quark condensate lj3S|l with u = 1. Line 4: Vanishing of four-quark condensate with 
u = 2. 



Note that here M denotes all mesons, including the Goldstone bosons. 

Now we can study the vanishing of the condensates also as functions of the energy and baryon density. Figure El 
shows the corresponding results. We observe first of all that the critical energy density widely differs for the four 
different lines. This is not surprising: Differences in the temperature multiply to a high power for the determination 
of the corresponding energy density (for massless states we would have e ~ T^). Figure O reveals a particularly 
interesting aspect: The energy density is extremely constant along the line where the four-quark condensate vanishes. 
Only for large baryon densities a different behavior sets in. This is the region of low temperatures (numerically 
T < 100 MeV) where we distrust the resonance gas approximation anyway. The critical energy density obtained from 
the vanishing of the four-quark condensate for u = 1 is about 0.8 GeV/fm'^. This agrees very well with lattice QCD 
results |37j . Note that the critical energy density obtained from the four-quark condensate for u = 2 is not quite that 
constant — albeit still more stable than the critical lines obtained from the two-quark condensate. We recall that 
u = 1 corresponds to our genuine estimate (|26|l whereas u = 2 has been introduced by hand to get an idea about the 
uncertainties. We conclude that the four-quark condensate seems to be a very useful order parameter to determine 
the line of chiral symmetry restoration — at least within the resonance gas approximation. 

Finallv we shall discuss the uncertainties induced by the fact that not all resonances listed by the particle data 
group Q are well established: We have checked that it makes no difference, whether we use all baryons of Ij or only 
the 3- and 4-star baryons. Figure 01 illustrates that it also makes no difference, whether one uses all mesons or only 
the lowest pseudoscalar and vector meson nonets. Using only the lowest octet and decuplet baryons instead of all 
baryons makes a difference, especially in an intermediate range of chemical potentials (see figure 0J. We conclude that 
high-mass baryons are to some extent important due to their large multiplicities whereas high-mass mesons are not 
relevant. We recall that this is a fortunate situation, since we do not know the overlap of high-mass mesons with the 
corresponding quark currents which appears in (|S^ . 
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FIG. 3: Vanishing of two- and four-quark condensates in the plane of energy and baryon density for different values of u. Same 
line code as in figure |21 
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FIG. 4; Vanishing of the four-quark condensate (for u = 1) for difi'erent sets of resonances. Full line: all mesons, all baryons; 
dotted line: lowest two meson nonets, all baryons; dashed line: all mesons, lowest two baryon multiplets. 
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VI. SUMMARY AND OUTLOOK 

The present work consists of two parts: First, a motivation why the four-quark condensate (|14|l should be regarded as 
a useful order parameter of chiral symmetry restoration. Second, a calculation within the resonance gas approximation 
to determine how this four-quark condensate drops as a function of temperature and baryo-chemical potential and in 
particular where it vanishes. As our favorite prediction for the line of chiral symmetry restoration we regard the full 
lines in figures Inland 01 This yields a temperature T^ w 174 MeV at vanishing chemical potential and an energy density 
tc ~ 0.8 GeV/fm"^ which stays remarkably constant as a function of baryon density or baryo-chemical potential. 

The corresponding line which one obtains from the vanishing of the two-quark instead of the four-quark condensate 
docs not show the same stability for the energy density Cc. This seems to be in contrast to the work ^j\ where also 
the in-medium change of the two-quark condensate has been determined within a resonance gas approximation. Here, 
we have to point out that our estimate for the resonance sigma terms H27(l adopted from [21| differs from the estimate 
utilized in "iofl . In the latter work the approximation [41| 

ax ~ Ml — (43) 

mx 

with 0.9 < A < \.2 has been used for all resonances, no matter whether these are baryons or mesons. Such an 
approach does not agree with large- A^c counting rules |2J|- For baryons the sigma term should scale with N^, whereas 
for mesons it should be constant. This requirement is satisfied by our approach: The factor 3 — iVg in H26|l and H27|l 
is actually Nc — Ng since a baryon consists of Nc quarks |4j|. In contrast, in the approximation H43|) the sigma terms 
of baryons scale with 1/Nc instead of Nc since baryon masses — which appear in the denominator of (|43|l — scale 
with Nc- Note that the numerator in H43|l cannot scale with Nc- Otherwise the scaling law for the meson resonances 
would be violated. Due to that deficiency we prefer our estimate (|27ll . 

Next we would like to discuss our approach in a more general context: In principle, we do not expect that any 
of our order parameters completely vanishes for high densities/temperatures. Only if chiral symmetry was an exact 
symmetry of QCD, the order parameters would exactly vanish in the chirally restored phase. In reality with finite 
current quark masses we expect a sizable drop of the order parameters at the phase transition or crossover point. After 
this transition point the order parameter presumably levels off. We have sketched the expected behavior in figure 
El Clearly the exact behavior of the order parameters near the transition point is hard to extract from low-density 
expansions or more general from an approach which uses only the degrees of freedom which are relevant below the 
transition point. Therefore we should not expect to obtain a completely accurate description of the transition region. 
Nonetheless, if one has a model which still works reasonably well in the region where the drastic change of the order 
parameter sets in, the estimated transition point should not deviate too much from the real one. Such a scenario is 
sketched in figure [SI 

For the evaluation of the in-medium behavior of the four-quark condensate we have used two crucial approximations: 
First, a resonance gas to model the in-medium system, and second, large- iVc arguments to evaluate the four-quark 
condensate for the non-Goldstone bosons. Concerning the latter approximation it is unfortunately hard to check its 
validity for the real world of A'c — 3. For general arguments in favor of the large- iVc approximation we refer to |43ll44l | 
and to the large amount of works which cite these seminal papers. 

Concerning the resonance gas approximation there is at least one aspect which might lead to additional non- 
negligible contributions beyond the ones obtained from the resonance gas: Since the pions are very light, the threshold 
for their interactions (2M^) is not significantly Boltzmann suppressed for the temperatures we are interested in (close 
to 200 MeV). On the other hand, if all the pion-pion interactions were mediated by narrow resonances, the relevant 
scale would not be set by the threshold but by the resonance masses. Indeed, there are arguments in favor of a 
resonance saturation of the low-energy constants which determine the pion-pion interaction (Sii. i49il iSOI]. However, 
these resonances, in particular the p-meson, are not narrow. Therefore a significant interaction already shows up at 
energies below the p-meson mass. Hence, pion-pion interactions should be taken into account beyond the resonance 
gas approximation. For the two-quark condensate at finite temperature and vanishing barro-chemical potential these 
interactions have been evaluated up to three loops within chiral perturbation theory in [Sj. It turned out that the 
effects induced by these interactions are not completely negligible, but small as compared to the combined effects from 
the non-interacting pions, nucleons and the resonances |3l|. To clarify whether this also is true for the four-quark 
condensate is devoted to future work. Another possible and straightforward extension of our approach is the inclusion 
of isospin and strangeness chemical potential. Also this is beyond the scope of the present work. 
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FIG. 5: Sketch of the real behavior of an order parameter as a function of temperature (full line) and of a reasonable approxi- 
mation obtained from the low-density degrees of freedom (dotted line). 
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APPENDIX A: RESONANCE DENSITY AND LARGE NUMBER OF COLORS 

From H29|) and (|30|l it seems that to leading order in \/Nc we only must consider baryonic states. However, there is 
an additional implicit iVc-dependence in px as we will discuss now. We will find that we have to distinguish different 
regimes for the baryo-chemical potential. This is intimately connected to the fact that the masses of baryons scale 
with Nr El: 



if X is a baryon, whereas 



mx = 0{N,) , 



mx = 0{N^) ., 



(Al) 



(A2) 



if AT is a meson. The reason is that a baryon consists of Nc quarks while a meson is a quark-antiquark state. 

We start with the case that the chemical potential is small compared to all baryon masses. For temperatures which 
are reasonable for hadronic degrees of freedom we have for baryons 



'mi^+k'^ - p> mx - ^^:$>T. (A3) 

In this case we can replace the Fermi distribution in (|23ll by a Boltzmann distribution and we can treat the problem 
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non-relativistically. This yields for baryons 

= o (iV3/2 exp(-iV,)) . (A4) 



For mesons there is no TVc-dependence, px — 0{N^). Since the expectation value of the four-quark condensate is 
enhanced by Nc for baryons (cf. (|29|l and HSOfl). we find in total that the relative importance of a baryonic contribution 
as compared to a mesonic one is o{Nc exp(— TVc)). In the large- iVc limit an exponential suppression overwhelms 
any power enhancement. Therefore, for small baryo-chemical potentials the in-medium change of the four-quark 
condensate is dominated by mesons. 

The situation changes, however, for larger baryo-chemical potentials to which we turn next: For mx — ^J^ — 0{N^) 
there is no implicit A'c-dependence left in px, i.e. px = 0{N^) even for baryons.^" In this case, the in-medium 
four-quark condensate is dominated by baryons. 

To get a serious estimate for all regions we should keep the respective leading order in l/Nc for every type of hadron. 
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